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Abstract. In this paper we are giving a quantitative description of two different configurations 
for noncontact gears. We consider the solutions from a perturbative calculation for two 
semitransparent parallel plates and concentric cylinders both with corrugations on the inner 
surfaces. In the case of corrugated parallel plates we discuss results from first- and second-order 
perturbation calculation in the corrugation amplitudes and we will concentrate on the first-order 
perturbation for the case of the corrugated concentric cylinders (the second order calculation 
is under study), both for the weak and strong couplings. We compare the perturbative results 
with the results from the PFA and an exact weak coupling calculation. 



1. Introduction. 

When two parallel plates with corrugations are situated close to each other, besides the usual 
Casimir attraction that takes place between them [1], there exists another contribution to the 
attraction that depends on the corrugations of the surface of the plates and therefore, on the 
non-normal coordinates (in the case of the parallel plates this coordinate is along the parallel 
surfaces where the corrugations are imposed, and in the case of the concentric cylinders this 
coordinate is the azimuthal angle with origin at the center of the cylinders). This new term 
gives rise to a lateral force or torque that could, in principle, make one plate slide with respect 
to the other or one cylinder rotate with respect to the other. Here we discuss such a situation. 
We base our discussion on a perturbative calculation details of which can be found in [2-4]. 
The case of QED parallel plates has been worked out in the perturbative regime in [5,6]. Some 
nanoscale mechanical devices based on the lateral Casimir force between different geometries 
have been proposed in [7-9]. 

We divide the paper in two relevant parts. The first one is devoted to corrugated parallel 
plates. We show and discuss the perturbative results for two limiting cases, strong and weak- 
coupling. The weak-coupling limit is also compared with an exact calculation and the PFA 
solution. The second part of the paper shows the results for corrugated concentric cylinders. 
We give perturbative results for both the Dirichlet and the weak-coupling limits. We extrapolate 



Figure 1. Parallel plates with sinusoidal corrugations. 



these results to the case of very large radii of the cylinders while keeping constant the distance 
between them. This situation corresponds to parallel plates. 

2. Parallel plates with corrugations. 

Let's consider two semitransparent plates lying on the x-y plane and placed parallel to each other. 
The plates have corrugations on the facing surfaces along the y-axes. We consider translational 
invariance on the x-axes and we work at a given frequency. The average distance between the 
parallel slabs is a and we model the corrugations by the functions hi{y) and h2{y). These are 
defined such that at any point the inequality a > hi{y) — /i2(y) holds. 

For the cases studied here we choose sinusoidal corrugations, hi{y) = hism{kiy). Each plate 
is described by a potential of the form, 

Vi{z,y) = Xi6{z - ai - hi{y)), (1) 

where i = 1,2 are labels that identify the individual plates. In this way we can consider our 
system as two semitransparent parallel plates making up the background and described by the 
potentials V^^\z) = Xi6{z — ai) for z = 1,2. On the top of it we have the corrugations hi{y) 
which can be considered to be a perturbation of the background (see figured]). 
Let's make one of the plates shift a distance yo with respect to the other one. 

hi{y) = /iisin[A;o(y + yo)], (2a) 
^2(y) = h2sm[koy], (26) 

where ko = 27r/d is the wavenumber corresponding to the corrugation wavelength d 0. This 
gives rise to a force in the lateral direction, the so-called lateral Casimir force, 

^ dE 

dyo 

Since the corrugations can be considered a perturbation over the background V^l'^\z), we can 
easily extract their contribution by calculating the deviation of the total potential Vi{z,y) from 
that of the background, 

AVi{z,y) = V,{z,y) - v/"' = \,S(z - ai - hfy)) - >.,S{z - ai) 

n=l n=l 



^ The same wavelength is assumed in both plates because it is only in this case that we find a leading order 
contribution to the lateral force. 



The key formula to calculate the Casimir energy is [10-12], 



AE= — TrlnGG(°) 
2r 



-1 



(5) 



where the Green's function G satisfies a differential equation with a potential function that 
describes the corrugations in both semitransparent plates, 



d'^ + Vi + V2 



G = 1, 



(6) 



and the corresponding Green's function associated with the background obeys the differential 
equation 



(0) 



V, 



(0) 



= 1, 



(7) 



where the corrugations are not present. However, since only the term containing the perturbative 
potential in both plates play a role in determining the lateral Gasimir force, we can reduce the 
above more general formula ([5]) to 



£^12 = Trln 

2r 



1 - G1AV1G2AV2 



(8) 



where now Gi refers to the Green's function when only one of the plates has corrugations on it. 
It obeys the equation 



(0) 



Gi 



1. 



(9) 



It is worth noticing that only the corrugations from plate i enter this equation. 

Expanding the above formula dS]) in a second- and fourth-order perturbative calculation in 
the corrugation amplitudes, gives expressions for the leading- and next-to- leading-order lateral 
Casimir energy. The general formulas turn out to be very long and complicated and they will 
not be shown here (details of the calculation as well as final formulas can be found in [2,4]); 
here we focus on the solutions for the Dirichlet and weak coupling limits. 

2.1. Dirichlet limit. 

Here we consider the strong coupling limit, where a\i^2 ^ 1) both in the leading and next-to- 
leading orders. 



2.1.1. Leading order. Expanding ([8]) to second order gives us the first non-zero contribution 
to the interaction energy. 



E 



(2) 
12 



LxLy 



(10) 



where the function A^^'^\kQa) has the form 

Jo 



dt- 



sinh s sinh Sj 



(11) 



and it is normalized such that ^^'^^(0) = 1, which becomes a convenience when we compare the 
results with those obtained by the proximity force approximation. We have used the notation 
= + t^, s\ = + {t ^ ^o)^) and to = ^o^J- This expression can be numerically evaluated 
and has been plotted in figure [21 
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Figure 2. Dirichlet limit: Plot of Ajj' {k^a) versus koa. 



The lateral Casimir force can be calculated using ([3]) which yields 
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(12) 



where l-^^calol magnitude of the normal scalar Casimir force between two uncorrugated 

parallel Dirichlet plates given as 



p(0) 
-^Cas,D 

LxLy 



' 480 ' 



(13) 



2.1.2. Next-to-leading order. The third-order term in the perturbation expansion gives zero 
contribution in the particular case when the wavelength is the same in both plates. If we do 
not restrict ourselves to this condition, the third order would give a finite contribution to the 
interacting energy becoming, in this case, the leading contribution itself. Since we are assuming 
ki = k2 = ko, we expand dS]) up to fourth order to get the next-to- leading-order contribution to 
Ei2. In the strong-coupling limit this gives us 



E- 



(4) 

12 



TT^ hi /l2 15 

240 a^~a~a^ 



cos(fcoyo)<! 4^§''^(^oa) + M Ag'=^)(fcoa) 



hi 

2 "^D 



cos(2A;oyo) \ — — A^'^^ (fcoa) 
1 a a 



(14) 



where we have introduced the functions, 
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tanhstanhs- tanh s tanh 
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tanh s sinh s+ sinh s_ 



(15a) 
(156) 



These functions can be calculated numerically and they are plotted in figure E] for a later 



because in such a case the second order expansion exactly cancels. 
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Figure 3. Dirichlet limit: Plot of A^^'^\kQa) and A^'^''(/coa) versus k^a. The dashed curve 
represents A^^'^\kQa) which is plotted here for reference. 



(2,2), 



discussion. 

The next-to-leading-order contribution to the lateral Casimir force in the Dirichlet limit reads 



-^Lat D = 2 ^oa sin(/coyo) 
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2.1.3. Discussion. In order to compare our results, we are going to confine ourselves to the 
case where the amplitudes of the corrugations in both plates are the same, /ii = /12 = h. Also, 
by observing equations (jl2p and (jl6p we notice that there is a common factor in all of them. It 
involves the normal Casimir force, as well as terms dependent on the geometry of the corrugations 
and on the initial shift 7/0 between the plates. By removing this common factor, what is left 
is just a dimensionless weight factor. We can therefore normalize the above expressions to the 
common factor and compare the resulting dimensionless term that we obtain in this manner 
from each of the mentioned formulas. In this spirit we define the weight factors as 



Td ( k^a, - ,koyo 
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(17) 



We shall use superscripts [m, n] on these dimensionless quantities to denote the order in the 
power series expansion to which the lateral force (not J-) has been calculated. In particular m 
will signify the order in the parameter koa, and n will denote the order in h/a. The perturbative 
results are complete in the koa dependence and thus will be denoted by jr[°°'"]. In these notations 
our results can be summarized as (see equations (|12p and (|16p ) 
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where means weight factor of the lateral Casimir force up to the fourth order. That is, it 

includes both leading and next-to-leading orders. 

We can in fact look at how much the next-to-leading order contributes to the force compared 
to the leading order. The left plots in figure [H show the second-order contribution to the force 



with a dashed hne and the fourth-order contribution with a solid hne for different values of the 
angle and constant values of koh. We see that for a fixed offset fcoyo the fourth order contribution 
becomes more important for high koh that means, for bigger corrugations. We show the plot 
for koijQ = to illustrate this fact more clearly because in this case the weight factor does not 
diverge for small k^a. However, at that value of the offset, the lateral force is zero as it can be 
noticed straightforwardly from ()12p and (1160 . For any other value of the offset, there are points 
where the graph diverges. This is because as koh approaches 2koa, the expansion is not longer 
valid. It corresponds to the case when the distance between the plates is smaller than the sum 
of the corrugations on both plates. 

On the right hand side of the same figure we plot the fractional correction of the next-to- 
leading-order contribution relative to the leading-order. It shows that at any given offset koyo 
and kQh the fractional error is higher for small separations and big corrugations give higher 
next-to- leading-order relative errors. 

In the absence of an exact answer, for the Dirichlet case, for comparison is not possible 
to extract the precise error in the perturbative results. In the next section we shall evaluate 
the corresponding result in the weak coupling limit and evaluate the lateral Casimir force non- 
perturbatively. We show that the error in the lateral force, for the weak coupling limit, after 
the next-to-leading order has been included is small for k^h <^ 1. Presuming that the same 
would hold for the Dirichlet case, we can expect the error to be sufficiently small to use the 
fourth-order results for comparison with experimental data. 



2.2. Weak coupling limit. 

We consider now aXi^2 to be very small and we study the second and fourth order in the 
perturbation series. 



2.2.1. Leading order. In this case the expansion of ([8|) leads to an expression that can be 
evaluated analytically 
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where we have introduced the function 



(0) 

Cas,W 



hi /l2 .(2)., X 

a a 



4(2)^7. N _ ''0 ^ 



tl 





— e 
to 



-to 



-to 



ml e*o 

m=0 



plotted in figure O We evaluate the lateral Casimir force to this perturbation order as 
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2.2.2. Next-to-leading order. We can write the total fourth-order contribution to the interaction 
energy in the presence of corrugations in the weak limit to be 
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where in contrast to the result in the Dirichlet limit we observe the appearance of the same 
function, though with different arguments, as coefficients. This function has been suitably 
defined as 



4! dt^ 



1 



-to 



-to 



E 

m=0 



64(^0 

pto 



(23) 






Figure 4. Dirichlet limit: The plots on the left show ' versus k^a for various values of 

k^UQ. The dashed curve represents ^]^'^^ which is plotted here for reference. In each plot the 
higher values of k^h deviate more from the reference. The plots on the right show the fractional 
correction in the next-to-leading-order contribution. In each plot the higher values of koh have 
larger corrections. 
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Figure 5. Weak coupling limit: Plot of A^{koa) = e ^""■e2{koa) versus kga 
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Figure 6. Weak coupling limit: Plot of (koa) and {2koa) versus koa. The dashed curve 

(2) 

represents A\y (koa) which is plotted here for reference. 



so that it equals unity at to = 0. These are plotted in figure [6l Using ([22]) in ^ we find that 
the fourth-order contribution to the lateral Casimir force in the weak coupling limit equals 
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2.2.3. Proximity force approximation. It is interesting to find this result since, as we shall see 
later, the weak-coupling limit can be solved exactly. That will allow us to compare both the 
perturbative and the proximity force approximation results to the exact solution. For sufficiently 
small distances a (in comparison to d = 27r /c^^) we can treat the plates to be built out of small 
sections for which the energy is approximately that of the parallel plate geometry. This is the 
so-called proximity force approximation. To compute the force under this approximation we 
consider the distance between the plates to be 



a{y) = a + /i2 sin[A;oy] - hi sin[A;o(y + yo)], 



(25) 



so that we can integrate the resulting energy after substituting the above in the expression for 
the Casimir energy between two plates in the weak coupling limit. 



Therefore, in the PFA, the lateral Casimir force in the weak coupHng Hmit reads, 
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2.2.4- Nonperturbative result. It was shown in [10] that exact results for the Casimir force, in 
the weak coupling hmit, can be achieved for specific geometries. This is being extended for a 
class of geometries, for the scalar case, in [13], and for dielectrics in electromagnetism [14]. In 
our case the solution of the exact weak coupling reads 
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2.2.5. Discussion. Plots analogous to the ones showed in figure U] can be drawn for the case of 
weak-coupling. However, since these have the same physical content as those for the Dirichlet 
case we do not show them here. Instead we are going to compare the perturbation expansion 
and the PFA to the exact weak-coupling solution. 

In order to do that, we follow the notation in (|18p and normalize the force for the weak- 
coupling as 
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The entire expression up to fourth order in the perturbation series in the weak coupling reads 
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The result in the proximity force approximation should be valid to the first order in koa and any 
order in h/a. Thus !F^^'°°^ will represent the PFA result in our notation. In the weak coupling 
this is 
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and for the exact calculation from ([27]) gives 
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The above non-perturbative expression for the lateral Casimir force on corrugated plates (in 
the weak-coupling limit), has been plotted in figure [71 The perturbative results are plotted as 
dashed curves in the same figure. We observe that the perturbative results, when the next-to- 
leading order is included, compares with the exact result remarkably well for koh <C 1 keeping 
the condition that 2/i < a. The fractional error in the lateral force when the next-to-leading 
order contribution is included is displayed on the right column in figure [71 For a fixed offset, 
the error gets larger for higher values of kph. For koh small, we get very good approximation in 
the whole range of koa. 

We can also compare the exact result with the one obtained from the PFA. In figure[8]we plot 






Figure 7. Weak coupling limit: The plots on the left show Tj-, ' versus k^a for various values 

of kono at k^h = 0.3. The dashed curves represents J^j^'^' and which is plotted here for 

reference. In each plot the exact value is greater than the value estimated by J^]^'^' . The plots 
on the right show the fractional error in the perturbative result. In each plot the higher values 
of k^h have larger errors. 
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Figure 8. Weak coupling limit: Fractional error in PFA is plotted as the dashed curve. The 
corresponding error in the perturbative result is plotted as the bold curve and is described by 
the axis on the right. 




Figure 9. Weak coupling limit: On the left the PFA is shown as the dashed curve. The exact 
results for higher values of koa are seen to deviate from the dashed curve. On the right the 
fractional error in PFA is shown to increase with higher values of koa. 



the fractional error in the PFA result, for fixed koh, and compare it to the fractional error in the 
perturbative results. We note that PFA is a good approximation for kQa <C 1 and perturbative 
results are valid for k^h <C 1. Both the PFA and the perturbative analysis are restricted to 
2h < a. 

We earlier noted that the PFA limit is obtained by taking the limit kf^a ^ while keeping 
the ratio h/a fixed. In figure [9] we plot the lateral force in the PFA limit and compare it with 
the exact result for various values of koa. We note that the error in the PFA is less than 
1% for koa <^ 1 for arbitrary h/a. For fcoyo = '^/'^^ we observe that the PFA is a very good 
approximation for h a which satisfies 2/isin(A;oyo/2) < a. After viewing the plots for various 
values of fcoyo we note that in general the PFA is a good approximation for /i ~ a and further 
beyond for offsets fcoyo < vr/4. It is, in fact, plausible that the PFA holds for large amplitude 



Figure 10. Non-contact gears: Concentric corrugated cylinders with the same corrugation 
frequency, v = 15, on each cyhnder. is the angular shift between the gears. 



corrugations for small offsets because the corrugations fit together like fingers in a glove. 

The experiments are highly grounded on calculations based on the proximity force 
approximation. There has been some misunderstanding regarding the range of validity of the 
PFA versus numerical calculations [6, 15, 16]. As we have shown here the ranges of validity of 
both solutions may not overlap. 

3. Concentric cylinders with corrugations. 

In this section we show the numerical calculation for the leading-order Casimir torque between 
two concentric corrugated cylinders. The cylinders have average radii ai and a2 and the mean 
distance between them is a = 02 — ai like it is shown in figure [TOj The set-up is equivalent to the 
one in the previous case. The cylinders are described by a potential of the same kind as before. 



where i = 1,2 identify the individual cylinders and hi{6) are the functions describing the 
corrugations. For an initial angular shift 9q as shown in figure [lOl and assuming sinusoidal 
corrugations, we have 



where hi^2 are the corrugation amplitudes and f is the frequency associated with the 
corrugations. 

In the next sections we will show the results for the strong and weak coupling limit up to 
second order perturbation. The Casimir torque is calculated by taking the derivative of the 
interaction energy with respect to the shifted angle. 



Details of the calculation can be found in reference [3,4]. In the last section we discuss the 
results. 

3.1. Dirichlet limit. 

In the Dirichlet limit the interaction energy can be expressed in the form 
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t igure 11. Dirichlet limit: Plots of bP^{^) versus to, for to < and fixed The dashed 
curve is the corresponding plot for corrugated plates which is approached by the corrugated 
cylinders for larger values of u. 



where we have divided by a factor of 2irR, which is the mean circumference. The function 
is defined as 
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where a = R= "^^"^ and the functions Dm are given by 

Dm{a.] x) = Im{x[l + a])Km{x[l - a]) - Im{x[l - Q\)Km{x[l + a]), 
Using ([M]) we evaluate the Casimir torque to be 
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3.2. Weak coupling limit. 

According to reference [3] the interaction energy in the limit a\i <C 1 becomes, 



where we have defined the function 
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The Casimir torque per unit area, for the weak coupling limit, can thus be evaluated to be 

r(2)^ _ _ A1A2 h h2 
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The function B^if^^ {a) is plotted in figure [T2l 




Figure 12. Weak coupling limit: Plots of BI {■^) versus to, for to < and fixed v. 
The dashed curve is the corresponding plot for corrugated plates which is approached by the 
corrugated cylinders for larger values of v. 



3.3. Discussion. 

We can compare our results with the results from the corrugated parallel plates. In the limiting 
case when the two cylinders have very large radii but the distance between them is kept constant 
we can consider that we have parallel plates in a region of small variations in the angle ^o- This 
corresponds to m, n — oo such that m/R\s finite and a — > 0. 

In figure [TT] we plot the function B^"^^ {a), where to = 2ai^. The dashed curve corresponds to 
the parallel plates. We see that the plot approaches the result for parallel plates as we increase 
V. Notice that a cannot get values bigger than 1 since that case corresponds to the inner radius 
equal to 0. When a — > we reach the PFA limit. 

The same conclusions can be drawn from figure [T^] for the weak coupling case. 

In general we see that our results for the Casimir torque reproduce the results for the lateral 
force on corrugated parallel plates in the limit of large radius and small corrugation wavelengths. 
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